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In many organic molecules the strong coupling of excess charges to vibrational modes leads to 
the formation of polarons, i.e., a localized state of a charge carrier and a molecular deformation. 
Incoherent hopping of polarons along the molecule is the dominant mechanism of transport at 
room temperature. We study the far-from-equilibrium situation where, due to the applied bias, 
the induced number of charge carriers on the molecule is high enough such that charge correlations 
become relevant. We develop a diagrammatic theory that exactly accounts for all many-particle 
correlations functions for incoherent transport through a finite system. We compute the transport 
properties of short sequences of DNA by expanding the diagrammatic theory up to second order in 
the hopping parameters. The correlations qualitatively modify the I-V characteristics as compared 
to those approaches where correlations are dealt with in a mean-field type approximation only. 



PACS numbers: 71.38.-k, 72.80.Le, 05.60.-k, 87.14.gk 

I. INTRODUCTION 

Molecular electronics experiments performed during 
recent years have probed the conductance and current- 
voltage characteristics of a large variety of molecules. 
Several experiments on long molecules indicate that 
transport is is not described by coherent Landauer trans- 
port or tunneling but rather by an incoherent hopping of 
charge carriers along the molecule. Examples are exper- 
iments on DNAii 2 - or oligophenyleneimine wires 3 -. In the 
latter experiment the length dependence of the conduc- 
tance clearly demonstrated the crossover from the coher- 
ent (tunneling) to the incoherent transport regime at a 
molecule length of about Anm. 

In many experiments, the molecule consists of repeated 
segments (either identical or with chemical modifica- 
tions) where the quantum mechanical hopping amplitude 
between the segments can be tuned to some degree by the 
choice of the "linker group". If the hopping amplitude 
between the segments is large (e.g., for stiff molecules 
with fully conjugated electron systems) the quantum- 
mechanical coherence on the molecule can extend quite 
far at low temperatures, such that even molecules of up 
to several nm length display signs of coherent transport, 
at least within the molecule 4 ^. On the other hand, if the 
coupling of segments is weak (e.g., for flexible molecules 
with weakly conjugated or saturated linker groups) the 
coherence decays quickly, such that charge carriers are 
typical localized over a single or a few segments only. 
In this case, the molecule tends to change its conforma- 
tion in order to lower its energy when charged, a process 
called polaron formation. The polaron is a combination 
of a charge carrier and the localized deformation. At 
room temperature, charge transport is then dominated 
by incoherent hopping of polarons along the molecule. 
At low temperature, coherent "band-like" transport of 
polarons might be observable. 

The theoretical description of polaronic effects in 



molecule-electrode setups so far either focused on molec- 
ular single-level systems^!^ or described polaron hop- 
ping in long molecules by assuming 'simple' rate equa- 
tion s ' n . Bottger and Bryksin have shown in a rig- 
orous description of polaron transport in bulk systems 
that even in the absence of Coulomb interactions phonon- 
mediated charge correlations between different sites de- 
velop^. However, in their calculations they included 
these correlations only in a mean-field like manner. The 
earlier rate equation approaches to polaron hopping in 
nanoscale systems also treat correlations within this 
mean-field approximation. 

The mean-field approximation to many-particle cor- 
relations can be justified in systems with low density 
of charge carriers, and is a sufficient approximation for 
many doped (organic) semiconductors. In molecular elec- 
tronics experiments, however, where a transport bias on 
the order of IV is applied over a molecule of a few 
nanometer length, the average charge density may be 
much higher, and correlations become relevant. For ex- 
ample, for small molecules often the Coulomb interac- 
tion (or charging energy) dominates, leading at low tem- 
peratures to transport characteristics similar to single- 
electron transistors ! 13 ! 14 With increasing molecule size 
the relevance of Coulomb blockade decreases, but still, 
the transport along the molecule is affected by charge cor- 
relations, either due to (non-local) Coulomb interaction 
or the (retarded) interactions mediated by the coupling 
of charge carriers to vibrational modes mentioned above. 
We will demonstrate that in general such correlations are 
not sufficiently described by a mean-field approach. 

We have extended the diagrammatic approach by 
Boettger and Bryksin to describe molecular systems cou- 
pled to metallic electrodes. In the usual diagrammatic 
approach to small molecules (or quantum dots) 15 ! 16 the 
molecular eigenstates including the Coulomb interaction 
are the basis of a perturbative expansion in the weak cou- 
pling to metallic electrodes. In contrast, in the present 
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problem the "basis states" are "local" to the molecule 
segments (e.g. a DNA base pair or a phenylene ring). 
The expansion parameters also include the small hop- 
ping amplitudes between the molecular segments. As the 
perturbation expansion usually involves a "self-energy" 
resummation, where diagrams of a certain type (but of 
arbitrarily high order) are accounted for, the occupation 
numbers (or, in general, one-particle correlation func- 
tions) are coupled to higher order correlations functions, 
unless the system is non-interacting. Similar to, e.g., 
equation-of-motion methods, a hierarchy of equations for 
the correlation functions can be generated, which has to 
be truncated in some manner (often following more nu- 
merical necessities than physical arguments). 

The present paper demonstrates that an exact descrip- 
tion of correlation effects mediated by vibrational modes 
is possible for a finite size system. This is because the 
hierarchy is naturally truncated at the level of the high- 
est possible correlation function involving all 'sites' of the 
molecule. The resulting finite set of coupled linear equa- 
tions for the occupation number and the many-particle 
correlation functions can then be solved without resort- 
ing to any further truncation procedure. 

As an example we study short DNA molecules coupled 
to metallic electrodes. The DNA is modeled by a tight- 
binding chain identifying each base pair with a single 
tight-binding site. The sites describing either guaninc- 
cytosine (GC) or adenine-thymine (AT) base pairs have 
different onsite energies and are coupled by direction- 
sequence-) dependent hopping integrals (compare Ta- 
ble P). The polarons are formed by strong coupling of the 
charge degrees of freedom to local base pair vibrations. 
To ensure energy dissipation, these base pair vibrations 
are in turn coupled to a set of harmonic oscillators, de- 
scribing the influence of a dissipative environment. 

In the first part of this article we introduce the dia- 
grammatic technique describing incoherent polaron hop- 
ping transport through molecules or other nanoscale sys- 
tems which are coupled to metallic electrodes. This tech- 
nique allows the description of polaron transport with the 
exact consideration of correlation effects arising from the 
electron- vibration interaction. The approximation of the 
technique lies in the need to restrict the order of the ex- 
pansion in small hopping parameters. In the second part 
we apply this diagrammatic technique to polaron hop- 
ping transport through short DNA molecules coupled to 
metallic electrodes with the following results: (i) Correla- 
tions effects beyond the mean field approximation are rel- 
evant for the linear conductance in inhomogeneous DNA 
molecules already for low charge densities of less than 
1 percent, (ii) When a transport bias is applied over 
the molecule, correlation effects become important even 
when the occupation in equilibrium is negligible, (iii) In- 
homogeneous DNA molecules in general exhibit two max- 
ima in the zero bias conductance as a function of equilib- 
rium chemical potential (gate voltage) and also in the dif- 
ferential conductance as a function of applied transport 
bias. In contrast, the mean field approach only displays 



one maximum, (iv) Depending on the sequence, the sec- 
ondary maxima can be suppressed, as a consequence of 
a small hopping rate limiting the transport through the 
system. Details of the diagrammatic technique are shown 
in the Appendices lAllCl 



II. MODEL AND TECHNIQUE 

The minimal Hamiltonian to describe polaron trans- 
port through DNA is H = H el + H vih + H cl _ vib + H L + 
Hr + H t ,l + H t ,r + -ffbath with 

i <ij> 

At,l + #t,r = ^ [^cjjr^ 4- t r a\c nr 

a i ^ ' 

ffel-vib = Y Xai ° l ( Bai + B m) ■ W 

a i 

The term H e \ models the electrons on the molecule 
with operators oj, aj in a single-orbital tight-binding rep- 
resentation. This implies that the molecule consists 
of N parts (labeled i). The electronic properties of 
the molecule can then be described by the molecular 
orbitals (usually the HOMO or LUMO) of these sub- 
entities with on-site energies ej and hopping tij between 
neighboring parts of the molecule. The terms -Hl/r 
refer to the left and right electrodes. They are mod- 
eled by non-interacting electrons, described by operators 
c 1l/r> c «l/R: with a flat density of states p e (wide band 
limit) . Since we do not focus on the details of the cou- 
pling between the molecule and the electrodes, it is suf- 
ficiently described by 7?t,l + -Ht.r- The tunneling am- 
plitudes are assumed to be independent of the molecu- 
lar orbital i and the quantum numbers of the electrode 
states v. The coupling strength is then characterized by 
the parameter T L ' R cx p e \t L > R \ 2 . 

Polarons are formed due to strong coupling of elec- 
tronic and vibrational degrees of freedom. The vibrations 
labeled a, i are described in ff v ib> with bosonic operators 
B a i and B ai for the vibrational mode with frequency 
uj a i, i- e. every part of the molecule can vibrate indepen- 
dently. i/ i_ vi b couples the electrons on the molecule to 
the vibrational modes, where X a i is the strengths for the 
local electron-vibration coupling for the site i and mode 
a, respectively. To ensure energy dissipation and ther- 
mal occupation of the vibrational states we couple every 
vibration to its own bath if^bath, the microscopic details 
of which do not matter. 

A perturbative treatment of the strong electron- 
vibration coupling in the above Hamiltonian is not rea- 
sonable. Nevertheless, to allow for a perturbation ex- 
pansion we apply the so-called polaron or Lang-Firsov 
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unitary transformation 

H = e s He- s 

with the generator 



a] en 



frhJn 



B a i — B 



(2) 



(3) 



We introduce transformed electron and vibrational oper- 
ators, 



a-iXi 



R R ~t„ 



huj n 



and polaron operators 



exp 



(4) 



Operators \i with different indices i act on different vi- 
brational states, therefore they commute for all times. In 
terms of these quantities the Hamiltonian reads 



H =H Q + H' 
Ho = ^2( e i~ \)a} a i + ^2 ^ }hJai ( B ai Bai 



H' 



A; 



+ H L + H R 

: - ^ a hUjXj 

<ij> 

+ ^2 [ trc ir a tXt + t r *a\x\cv 
X 2 



hhJ 



(5) 
(6) 



physics in this regime we extend a formalism developed 
by Bottger and BryksirJ^ for polaron transport in bulk 
systems to account for coupling to metallic electrodes. 
To calculate quantities of interest, e. g. the occupation 

number l^a\{t)ai(t)^ and the current in a non-equilibrium 

situation with applied bias, we make a real time expan- 
sion of the occupation number along the Keldysh con- 
tour. The evolution in the interaction picture introduces 
the time dependence 



Oi(t) = 

Bi(t) =B t e 



-i{ti-Ai)t 



ate 



-ieit 



-iuJit 



From here on we will use the shifted onsite energy = 
e, — Aj in all expressions. 

The occupation number of the molecule can be written 

as pi(t) = (^aj(t)ai(t)^ _ . We express it in the interaction 

picture, assuming that the perturbation H' is adiabati- 
cally turned on from the time to = — oo, 

Pi(t) = (ut (t,-™)alaiU 6o {t,-<x>))_ 

\ o ' Ho 

with time evolution operator 



U Ho (t, -oo) = T ^ exp 



dtH'r, (t) 



(7) 



A Taylor expansion of the time evolution operators 
in H' defines a diagrammatic expansion. The forward 
time-evolution operator U Ho {t,—oo) is expanded on the 
upper branch of the Keldysh contour, whereas the back- 
ward time-evolution operator (t, — oo) is expanded on 

Ho 

the lower branch (see Fig. [T]). The index Hq indicates 
that these operators are written in the interaction pic- 
ture. The time ordering operator 'T' in Eq. [7] (anti-time 



H' describes the perturbation to the exactly solvable 
Hamiltonian Hq. The perturbation consists of the hop- 
ping along the molecule and to and from the electrodes, 
where the operators \i account for the creation and ab- 
sorption of vibrations in the hopping processes. The per- 
turbation can be considered small if either the hopping 
strengths tij or t r are small and/or if the polaron binding 
energy A is large, as the terms in a perturbative expan- 
sion are proportional to 1/A. 



— OO 



-oo 



Ufj {t, -oo) 



U Ho {t, -oo) 



A. Real-time expansion 

There are two limits to polaron transport, coherent 
band-like transport and incoherent hopping transport. 
For weak electron-vibration coupling and low tempera- 
tures coherent transport dominates, whereas for strong 
coupling and high temperatures transport is a sequence 
of incoherent hopping processes. In this work we will 
focus on incoherent polaron hopping. To describe the 



FIG. 1: Schematic drawing of the Keldysh contour and the 
forward and backward time-evolution operators. The open 
and crossed circle (the clamp) represent the two operators O; 
and a], respectively, which are evaluated at time t. 

ordering operator 'T') ensures that the different times ti, 
arising from the Taylor expansion of the forward (back- 
ward) time evolution operator, are ordered in the correct 
way along the contour. Note, oftentimes forward and 
backward time evolution operators are combined and a 



contour ordering operator 'TV is introduced to ensure the 
correct ordering of times along the Keldysh contour i 17 i 18 

In performing the expansion in the time evolution op- 
erators, we obtain certain operator products, which we 
have to average thermally. Since Ho is quadratic in the 
fermion operators, these can be treated using Wick's the- 
orem. On the other hand, the vibrational operator prod- 
ucts, involving various operators Xi(fj)t cannot be fac- 
torized. The rules for the evaluation of these operator 
products are given in Appendix [C] 

A specific term in the Taylor expansion, is represented 
by a diagram with a certain number of vertices on the 
upper and lower branch of the Keldysh contour, where 
each vertex is proportional either to tij (a hopping ver- 
tex) or V (a tunneling vertex). Each vertex consists of 
one open circle O (symbolizing a destruction operator) 
and one crossed circle (J) (symbolizing a creation oper- 
ator). All circles belonging to operators acting on the 
molecule are drawn on the inside of the contour, whereas 
circles belonging to electrode operators are drawn on the 
outside of the contour (compare e.g. Fig. [2J. The dif- 
ferent vertices are connected by fermion (solid) and vi- 
brational (dashed) lines and belong to different times ti, 
which have to be (anti-) time ordered along the (lower) 
upper branch of the contour. 

A feature of this expansion is that certain diagrams 
are diverging even in first order. These diagrams can 
be identified by so called "free sections" (indicated by 
the dotted lines in Fig. [2]). A free section is a part of 
the diagram between two vertices (except for the clamp) 
where a vertical line can be drawn such that only inter- 
nal fermion lines are cut (neither a phonon line or an 
external fermion line belonging to the electrodes). In 
such a case, the vertical line always cuts an even num- 
ber of internal fermion lines, as many left- as right-going. 
These left and right-going lines are pairwise associated 
with the same site. In the evaluation of such a diagram 
this leads to a divergence^ Thus an infinite number of 
diagrams has to be summed up in a way similar to a 
'ladder'-approximatio n 12 i 20 i 21 . The regions in between 
free sections (excluding the clamp) are called irreducible 
blocks. They do not diverge. 

In Figure [5] such a ladder-summation of a second or- 
der irreducible block representing a tunneling process is 
shown. Similar to a Dyson series, the summation of an 
infinite number of diagrams can be written as a self con- 
sistent equation for the occupation number p. 

For the interacting system we consider, certain dia- 
grams lead to equations coupling the occupation number 
to many-particle correlation functions. This is illustrated 
in Figure [3] The dotted vertical lines denoted (i) and 
(ii) indicate free sections which lead to divergences. Let 
us concentrate on the free section (i) in Fig. 3(a) To 
cure the divergence due to (i) a ladder summation has 
to be performed over all diagrams with the same diver- 
gence, i. e. the same free section. This can be done in 
a complete and tractable manner by the introduction of 
many-particle correlation functions. 
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FIG. 2: Ladder summation of a second order irreducible block 
representing two tunneling processes. The irreducible blocks 
are separated by free sections (indicated by dotted lines) . The 
summation can be recast into a self-consistency equation for 
the occupation number p. 



Aside of this technical argument, there are also sim- 
ple physical arguments why (and when) the two-particle 
correlation functions affect the behavior of the occupa- 
tion numbers (and the current). Consider a hopping 
process from sites m to site n: the hopping probabil- 
ity is determined by a second order irreducible diagram 
and the occupation of the initial and final site. Two- 
particle correlation functions express the probability to 
find e.g. the initial site m occupied and the final site n 
empty, such that the hopping process can succeed. In 
general, the occupation of different sites is correlated, 
except in the trivial case when there is exactly one par- 
ticle on the (central) system. If the charge density is 
finite, but very low, the charges are well described by 
a Boltzmann distributio n 22 ^ 3 and the two-particle cor- 
relation functions can be factorized in a "Hartree-Fock" 
type of approximation. This approach was taken in ear- 
lier works i 11 ' 21 ! 2 ^ 3 - 



B. A hierarchy of many- particle correlation 
functions 



By inspection of the part of diagram Fig. 3(a) 
left of the free section (i) one notices that this 
resembles a diagram arising from the real-time 
expansion of the two-particle correlation function 

(ut^ajarnaUaiUfj^)), see Fig. [3(b)} Straightfor- 
ward generalization shows that the ladder summation for 
the diagram 3(a) consists of all diagrams that arise from 



the real-time expansion of this particular two-particle 
correlation function, placed to the left of the free sec- 
tion (i) , as indicated by Fig. 2] In this way, an infinite 
number of diagrams to the occupation number p can be 
accounted for by involving this particular two-particle 
correlation function with the particular irreducible block, 
given by the right part of the diagram Fig. 3(a) from 



the free section indicated by (i) to the last vertex before 
the clamp. Other irreducible blocks involve other type 
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FIG. 3: (a) Sixth order diagram in the real-time expansion of 
the occupation number (iJg (t)ajaiU # (t)^ (b) Fourth order 
diagram in the real-time expansion of the two-particle correla- 
tion function (Ufj {t2)aj fflmim a i^H (^ 2 ) V F ree sections are 
indicated by the vertical dotted lines. 



of many-particle correlations functions. What kind of 
correlation function is needed can be read off from the 
vertices of the irreducible block (following the contour) 
that are connected to the fermion lines crossing the free 
section (see also App. |B|) . 



Fig. |3(b)| is only a particular diagram to the 
expansion of the two-particle correlation function 

{t%) a \ a rnQm a lU H {t2)} ■ However, the arguments 
for replacing the free section (ii) in this diagram work the 
same way as for the expansion for the occupation num- 
ber p. Therefore, we can write a linear equation for this 
two-particle correlation function involving all other two- 
particle correlation functions, the occupation number p 
and other many-particle correlations functions. These 
many-particle correlations fulfill yet another set of linear 
equations. 



In this manner a hierarchy of linear equations is es- 
tablished. In principle, the approximation to the exact 
solution of the problem lies so far solely in the limited 
number of irreducible blocks that can be considered in a 
real calculation, i.e. in the order of expansion of the irre- 
ducible blocks in the hopping and and tunneling vertices. 
In practice, the question arises whether the hierarchy of 
equations can be solved exactly, or whether other ap- 
proximations need to be applied, like a truncation or a 
factorization of the many-particle correlations functions. 




FIG. 4: A second order hopping diagram. The full lines rep- 
resent fermion lines. The dashed line represents the sum of 
all possible vibrational lines arising from the diagrammatic 
rules, here it has a value F t + (ti — t2)F+ 2 {t\ — £2) + F^(ti — 
t2)+F+ 2 {t 1 -t 2 ). 



C. Truncation of the hierarchy 



An important feature of the real-time expansion of 
the correlation functions is that certain diagrams van- 
ish due to the Pauli exclusion principle. For exam- 
ple, the diagram depicted in Fig. [5] arises in the ex- 
pansion of the two-particle correlation function and de- 
scribes hopping process from site m = ni2 to 7713 at 
time ti and back to m at time t\ . To its left the ir- 
reducible block is coupled to the three particle corre- 



lation function 



,t 



with mi = I. For = mi = I there arises a special 
situation, the three particle correlation function is zero, 
since for fermions hi(l — nj) = with hi = a,ai. 



Furthermore, by similar arguments the Pauli exclusion 
principle leads to a natural truncation of the hierarchy of 
equations for any finite system. The A^-particle correla- 
tion function (u^ (Q)a\aia\a2 ■ ■ ■ ffljyGjvi7g- (0)^ can not 
couple to any higher order correlation function because 
they all vanish (recall that N is the system size.) Thus, 
a closed set of linear equations for the real-time expan- 
sion of all correlation functions can be constructed. The 
formal solution of this set of equations is simply a matter 
of back-substitution. 

In the earlier work a 11 ! 21 ! 22 ' 23 , a Hartree-Fock factor- 
ization was applied in which products of particle num- 
ber operators are replaced their expectation values. This 
leads to terms in the expansion that should not exist if 
the exact correlation functions were considered. The dif- 
ferences to the full theory are small, if the electron or 
hole densities are so low that they can be described by 
Boltzmann statistics. However, there is the additional 
complication that the factorization leads to a non-linear 
self-consistency equation for the occupation numbers pi. 
At finite bias it can be quite difficult to find a converg- 
ing solution, especially as the system size becomes larger. 
This difficulty is avoided in our present theory where a 
solution to the linear equation set can be readily found 
by standard numerical methods. 
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FIG. 5: A second order hopping diagram. The full lines rep- 
resent fermions, the dashed line represents the sum of all pos- 
sible vibrational lines arising from the diagrammatic rules, 
here it has a value F+(t\ — t2)F+ 3 (t\ — £2) + F^(ti — t^) + 
F+(ti-t 2 ). 



Summarizing the above, the relation 



dti 



E 



^fel^JWmilftl ~t) 



,{ti)W mima i(ti -t) + 



(8) 



for the time derivative of the occupation number is ob- 
tained. The diagrammatic rules for construction and 
evaluation of irreducible blocks W are listed in ap- 
pendix [X] The generalized (one- and two-particle) 
correlation functions g mi (g mi m 2 ) represent any or- 
der of the creation and destruction operators a mi , 
(a mi , , a m2 , aj„ 2 ) that arises in the free sections. Us- 
ing the commutation relations for fermions, all general- 
ized correlation functions g mi ...m- can be expressed by 
a sum of correlation functions p m - L ...m i (of the same or 
lower order) where we fix the order of creation and de- 
struction operators such that the creation operators at a 
site are to the left of there destruction counterparts, i.e. 
Pm 1 ...m i = ttJnjdm! . . .a* m .a mi . In the rest of this article 
we only use these ordered correlation functions p mi ... mi 
(note that the occupation number is naturally denned by 
the one-particle correlation function p mi )- 



D. Explicit equations for the considered model 

In this article we want to calculate the stationary state 
of the system when a constant, finite bias is applied be- 
tween the metallic electrodes. Therefore, all correlation 
functions must be constant in time as well, such that 
Eq. [5] reduces to 

d f* 

dtiW mx m 2 l(tx) + . . . . 

mi m 2 J 

(9) 

Furthermore, we restrict ourselves to the lowest non- 
vanishing order of diagrams, i. e. second order. Examples 
of two such diagrams in the expansion of the occupation 



number and two particle correlation functions are given 
in Fig. O and 2] When comparing these two diagrams it 
is clear that they are mostly identical except for two ad- 
ditional fermion lines in Fig. [5] which enter as Kronecker- 
delta symbol. Of course, the corresponding irreducible 
blocks W couple to different correlation functions deter- 
mined by the fermion lines leaving to and entering from 
the left. 

According to the 'mirror' rule, to every diagram one 
can construct its complex conjugate by moving a vertex 
from the upper part of the contour to the lower one (and 
vice versa) . Following the rules given in the appendix the 
two diagrams (Fig.[5]and[4]) together with their respective 
complex conjugates and including the time integral of 
Eq. [5] have the values 



m 1 m 2 1 



fa* 



tl 2 



dt e h 



x [Fi + (t)F+ 2 (t) + F+(t) + F+ 2 (t)} S lmi 



(10) 



and 



VV. 



x [F+(t)F+ 3 (t) + F+(t) + F+ 3 (t)] S lmi 5 mm2 . 

(11) 

Explicit expressions for the constants /Q and the func- 
tions (t) are given in the appendix [X] The common 
element of all second order hopping diagrams in the ex- 
pansion of the various correlation functions may be de- 
fined as 



x [F+(t)F+(t)+F+(t)+F+(t)] . (12) 

For the second order diagrams describing the hopping to 
and from the electrodes there exist also such common 
elements. They are respectively 



w 



wx 



2nh 
dE 
2^h 



f L (E))Cl (F+(E-h)+5{E-h)) 



(13) 



where T L = 2-K\t L \ 2 p e , fh(E) is the Fermi function in left 
lead, and F^~(E) is the Fourier transform of F^~(t). For 
the right interface a similar expression holds involving 
./,;: /•. : and F+(E). 

Evaluating all second order diagrams simplified rate 
equations for the single-particle occupation number can 
be stated 

^Pl = [ - (Pi - Plm) Wlm + (Prn ~ Pml) W mi (14) 



dt 



dt Pl 



{Pl - P12) VV12 + (ft2 - p2l) W21 ■ 



(15) 
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For the two-particle correlation functions on the inside of 
the molecule and connected to the left electrode the rate 
equations have similar form 



5'-XY-3'(all in eV) 



Jt Plr ' 



E 



(Plm ~ Plmn) Win + (p mn Plmn 



(Plm - Plmn) W TO „ + (/)/ )W n 



and 



dt 



P12 



(16) 

- Pl2 Wh + (p 2 - P12) W\ 

- {Pl2 - P123) W23 + (Pl3 ~ P123) W32 ■ (17) 



All other rate equations not presented above are con- 
structed in an analogous way. This leads to a closed set 
of linear equations for all correlation functions up to the 
order of the system size N. 

The theory is current conserving as Eq. [T3] and [15] are 
the continuity equations for the occupation of site I and 
1, respectively, which equal zero in the steady state sit- 
uation we consider. Thus the current can be computed 
at any point of the entire system. For convenience, we 
choose to compute the current through the left lead given 
by 



h = e 



PiW^ + (l-pi)W_ 



(18) 



Note that all many-particle correlation functions drop 
out of this expression. They influence the current only 
via their effect on the occupation at the first site, p\. 



E. Polaron hopping transport in DNA 

In this section we apply the above presented diagram- 
matic approach to polaron hopping transport in DNA, 
which we have already studied earlier in a more simpli- 
fied approach, involving the Hartree-Fock factorization. 11 
DNA consists of a sequence of base pairs adenine and 
thymine and guanine and cytosine, which form a dou- 
ble helical ladder structure. The electronic properties of 
DNA are determined by the HOMO (highest occupied 
molecular orbital), situated on the guanine and adenine 
bases, and the LUMO (lowest unoccupied molecular or- 
bital), situated on cytosine and thymine basest 

To describe polaron-hole hopping the relevant molecu- 
lar orbital is the HOMO. We describe the DNA chain in 
a minimal tight-binding model where each tight-binding 
site represents one HOMO either on a guanine or an ade- 
nine base. Both on-site energies ti and hopping integrals 
tij depend on the base pair sequence. For the direction- 
dependent hopping matrix elements iy we use the values 
obtained from density functional theory by Siebbeles et 
aii2£ Adapting these values to our model of base pairs 
we obtain the next-neighbor hopping elements listed in 
table |H 



X\ Y 


G 


C 


A 


T 


G 


0.119 


0.046 


-0.186 


-0.048 


C 


-0.075 


0.119 


-0.037 


-0.013 


A 


-0.013 


-0.048 


-0.038 


0.122 


T 


-0.037 


-0.186 


0.148 


-0.038 



TABLE I: Hopping integrals tij taken from Ref. |25j and 
adapted to our model. The notation 5'-XY-3' indicates the 
direction along the DNA strand (see, e.g., Fig. lb in Ref. I24I.) 



As shown by Alexandre et alr 2 ^ and also other 
author a 27 ! 28 ' 29 polarons are formed on DNA molecules, 
although the size of these polarons is still controversial. 
Fits to the temperature dependence of the linear conduc- 
tance of experiments on long (> 1000 base pairs) DNA 
segments like Ref. Q] support the idea of small polaron 
formation with a local DNA distortioni 30 i 31 In this work 
we assume that the size of polarons is restricted to sin- 
gle DNA base pairs. Such small polarons are formed due 
to strong coupling of the electronic degrees of freedom 
to local vibrational modes of the DNA base pairs. Ex- 
emplary, we consider only a single vibrational mode per 
base-pair, the so-called stretch modes with frequencies 
hijji — 16meV for a GC base pair and huii — 11 meV for 
an AT base pahi^ 2 - 

The electron-vibration coupling strengths are chosen 
in such a way that the reorganization energy or polaron 
shifts (compare Eq.©, A A = 0.18 eV and A G = 0.47 eV, 
fit the values extracted from experiments and listed by 
Olofsson et al££. These values probably underestimate 
the effect of the solvent on the reorganization energy. 

To ensure energy dissipation and thermal occupation 
of the vibrations each base pair i is coupled to a local 
environment, -ff^bath, the microscopic details of which 
do not matter. This coupling changes the vibrations' 
spectra from discrete modes uji to continuous spectra, 



Di(u 



1 

7T 



»7<M 



(u - uji) 2 + rn{uj) 2 (uj + toi) 2 + ^(lu) 2 

I 



(19) 



with frequency dependent broadening r]i(uj)& The actual 
form of rji(u)) depends on the properties of the bath. 
A reasonable choice which assures also convergence at 
low and high frequencies is r)i(cj) = rjo ^ 9(uj c — to) with 

7/0 = 0.5 meV and a cutoff of the order of huj c = 0.045 eV. 
To account for the spectral function of the base pair vi- 
brations, the substitution ^ Q — » J dioDi{ijf) has to be 
made in all equations introduced above. 



III. RESULTS 

The main focus of this article is the investigation of 
the influence of correlations on the electronic transport 



8 



characteristics. To do so, we compare our present theory 
with exact correlation effects to our previous approach in- 
volving a Hartree-Fock factorization that deals with cor- 
relations only on the mean-field level. For simplicity, we 
denote the latter approach as "mean-field correlations" . 

It should be noted that correlation effects do not influ- 
ence the transport properties of homogeneous DNA se- 
quences or other homogeneous molecules. This is because 
for a homogeneous system, not only pi m — p m i (this is al- 
ways true by definition of the correlation functions), but 
also the hopping rates for forward and backward hopping 
processes are identical, W; m = W m ;. Therefore, the two- 
particle correlation functions drop out of Equations [T4l 
and IT51 for the occupation numbers, and consequently do 
not influence the current. The I-V curves are therefore 
identical, no matter whether exact correlations or mean- 
field correlations are considered. 



A. Zero bias conductance 

As noted earlier the difference between the two ap- 
proaches is small, when the occupation of holes or elec- 
trons is small, so that the occupation is described by 
a Boltzmann distribution function. To demonstrate this 
Fig.[6]shows the zero bias conductance Go as a function of 
chemical potential fi for a DNA molecule with sequence 
AAAGAAAA with mean-field correlations (black solid 
line) and exact correlations (red dashed line). We have 
chosen of /i = to lie above the (polaron shifted) gua- 
nine and adenine states, i. e. in the HOMO-LUMO gap 
of DNA. As can be seen from the Figure, especially from 
the inset on logarithmic scale, for values of \i < — 0.5eV 
and above (j. > — 0.15eV the two curves agree with each 
other. In the first region the electron occupation number 
is very small, whereas in the latter the hole occupation 
is very low. On the other hand, in the central region 
of the plot both curves differ strongly. The zero bias 
conductance is much lower when correlation effects are 
exactly accounted for. Furthermore, the red curve ex- 
hibits two maxima around the chemical potentials that 
agree with the onsite energies of adenine (?a) and gua- 
nine (ec). The black curve, with mean-field correlation 
effects, only shows a single broad maximum. Similar 
to coherent quantum transport through molecules, the 
mean-field type approximation of correlations overesti- 
mates the (zero bias) conductance and current and only 
shows a very simplified energetical structure. 

As an illustration, the effect of exact correlations on 
the occupation of the bases for the sequence AAA- 
GAAAA is shown in Fig. [7] The curves show the dif- 
ference of the occupation for the two approaches for the 
first (black solid line) fourth (red dashed line) and eighth 
(blue dash-dotted line) base of the DNA molecule AAA- 
GAAAA as a function of chemical potential /i. The 
values were calculated for a small transport bias Vb = 
0.01 V, which is in the linear regime. The sign and mag- 
nitude of the occupation difference changes with the bias 
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FIG. 6: Zero bias conductance Go as a function of chemi- 
cal potential /i for the DNA molecule with sequence AAA- 
GAAAA with mean-field correlations (black solid line) and 
exact correlations (red dashed line). The parameters used 
are ea = — 0.26eV, eG = +0.25 eV relative to the zero point 
of the chemical potential, polaron shifts Aa = 0.18 eV and 
Ag = 0.47 eV, symmetric coupling to leads with linewidths 
Fl = Tr = 0.001 eV, vibrational energies Kuja = llmeV, 
huG = 16meV, and room temperature feT = 25meV. The 
inset shows the same plot on logarithmic scale to stress that 
both models agree for very low and high occupation. 




FIG. 7: Lower graph: Difference of the occupation for ex- 
act and mean-field correlations for the first (black solid line) 
fourth (red dashed line) and eighth (blue dash-dotted line) 
base of the DNA molecule AAAGAAAA as a function of 
chemical potential fi. The used transport bias is Vb = 0.01 V, 
which is in the linear regime. As reference in the upper graph 
the plot of Fig. [6] is repeated. Parameters as in Fig. [6] 



direction and strength, respectively. 

The maximum in occupation difference for the adenine 
and guanine bases is reached for chemical potential values 
which are close to the onsite energies of either adenine or 
guanine, respectively. The effect of the correlation is the 
strongest for the isolated guanine base in the center of the 
sequence, where the maximum difference in occupation is 
three orders of magnitude greater, than for the adenine 
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-1.5 -1 -0.5 0.5 1 1.5 

V b (V) 

FIG. 8: Differential conductance dl/dVb (logarithmic scale) 
as a function of applied bias Vb for the DNA molecule with se- 
quence AAAGAAAA, for mean-field correlations (black solid 
line) and exact correlations (red dashed line). Parameters as 
in Fig. [6] but with the equilibrium chemical potential of the 
electrodes fixed at n — OeV. 

bases. The occupations of the other adenine bases (not 
shown), are similar to the ones depicted in Fig. [71 



B. Finite bias differential conductance 

In experiments a variation of the chemical potential is 
rather difficult to achieve. In molecular electronics this 
is usually achieved by using a back gate electrode under 
an insulating substrate. However, for DNA the complica- 
tion arises that its conformational structure is much in- 
fluenced by the surface potential of the substrate and the 
electric potentials of the back gate. Oftentimes, though, 
the I-V characteristics can be probed in setups where 
the molecules are in free suspension (mechanical break 
junctions) or standing upright in molecular monolayers. 

When applying a transport bias over the molecule, the 
occupation of the various molecular segments (base pairs 
in the case of the DNA) depends on the magnitude of 
the applied bias "felt" at the particular location. The 
corresponding potential profile can be interpreted as a 
"local chemical potential" that changes with the applied 
bias. As this local chemical potential moves over the en- 
ergies of the base pair levels, similar effects as in Fig. [5] 
are expected. In Fig. [8] we show the differential conduc- 
tance dl/dVb as a function of applied transport bias 14,. 
The red curve including exact correlation effects shows 
two maxima for both positive and negative bias, whereas 
the black line with mean-field correlations only has sin- 
gle peaks. For small bias both curves agree very well, 
as the charge carrier occupation in this regime is very 
low. Note that again the mean-field approach overesti- 
mates the (differential) conductance everywhere, except 
for large positive bias, where the second peak exists in a 
bias region where the mean-field approach shows expo- 
nentially small conductance. 



From these calculations it becomes clear that at finite 
transport bias the charge carrier occupations have value 
regions, where mean-field type of correlations are insuf- 
ficient to describe polaron hopping transport through 
molecules like DNA. 

We now discuss the position of the maxima in the zero 
bias and differential conductance. For a homogeneous 
molecule, both Go and dl/dVh only have a single maxi- 
mum, namely when the chemical potential (or the trans- 
port bias) is in resonance with the level energy of the 
DNA base pairs (e G = -0.22 eV and e A = -0.44 eV). 
Including exact correlations, for the DNA molecule AAA- 
GAAAA there are two maxima in the zero bias con- 
ductance as a function of chemical potential, or in the 
differential conductance, for both positive and negative 
transport bias (see Figure [5] and [5} . These maxima can 
also be associated with the level energies of the two dif- 
ferent types of bases, adenine and guanine. However, 
the exact positions of the maxima in the zero bias con- 
ductance deviate slightly from their expected resonance 
positions due to charge rearrangement effects between 
guanine and adenine bases. Since the hopping rates in- 
crease strongly with temperature, the charge rearrange- 
ment also increases. Thus the position of the maxima is 
temperature dependent. 

For the differential conductance the positions of the 
maxima are shifted more strongly as a finite transport 
bias will lead to much stronger charge rearrangements 
("polarization"). This effect is increasingly important 
for the maxima at higher bias, i. e. the second maxima 
are shifted more strongly from the anticipated resonance 
positions of the adenine bases energy as compared to the 
first maxima relating to the guanine energies. The charge 
rearrangement and therefore the position of the conduc- 
tance maxima is quite sensitive to the considered DNA 
sequence (see also discussion in Ref. [TTh . 



C. Sequence effects 

Are there always as many maxima as different species 
of molecular segments with different onsite energies, if 
exact correlations are considered? The black solid curve 
in Figure [5] shows the zero bias conductance as a func- 
tion of chemical potential for the DNA sequence AAAG- 
GAAA. This sequence is only a slight modification as 
compared to AAAGAAAA studied above, nevertheless 
the zero bias conductance exhibits only a single maxi- 
mum, the 'adenine '-peak is missing. A similar behavior 
is observed for the differential conductance as a function 
of applied bias. 

Does this mean that for this sequence the adenine res- 
onance does not occur? It turns out that it is just sup- 
pressed. The red dashed curve in Fig. [5] is calculated for 
the same sequence, but with a hopping rate Wgg en- 
hanced by a factor 100. In this curve again two maxima 
are visible, both in the zero bias and differential conduc- 
tance. Obviously, the adenine resonances were present 
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AAAGGAAA 
AAAGGAAA; W *100 




FIG. 9: The zero bias conductance Go as a function of chem- 
ical potential p for the DNA sequence AAAGGAAA (black 
solid line) with exact correlations. Only a single maximum 
around the guanine energy is displayed. For the dashed line 
the hopping rate between the two central guanine bases is en- 
hanced by a factor 100. This shows that the maximum at the 
adenine energy in the sequence AAAGGAAA is suppressed 
due to the peculiar hopping parameters relevant for DNA. 
Other parameters as in Fig. [6] 



but the resulting transport is strongly suppressed by the 
low hopping rate Wgg for our model of DNA parameters. 

At first glance, it might appear strange that the hop- 
ping rate between two guanine base pairs is so much lim- 
iting the transport at energies related to the adenine base 
pairs. However, this effect can be easily explained by 
looking at the current between the two guanine bases of 
the above DNA sequence. From Eg . [X4l one can deduce 

I =e(p4- p i5 ) W45 + e (pg - P45) VV54 
= e(p 4 -P5)W G G , (20) 

where the sites 4 and 5 denote the two guanine bases and 
for the second line W45 = W54 = Wgg wa s used. 

As the current is conserved, the above equation is iden- 
tical to the current obtained from Eq. [18] Obviously, a 
small rate Wgg leads in general to a small current. How- 
ever, the important matter lies in the difference of occu- 
pations (p4 — ps) of the two guanine bases (for a small 
applied transport bias) that varies strongly depending on 
the chemical potential. For the DNA sequence AAAG- 
GAAA the current reaches its maximum value when the 
difference in occupation between the two guanine sites 
is the greatest. This is the case when the guanine on- 
site energy is in resonance with the chemical potential, 
i.e. around p = — 0.22eV, as at this energy the guanine 
occupation changes rapidly from unity to zero while low- 
ering the chemical potential. For much lower values of 
the chemical potential (in particular p w — 0.44eV) the 
occupation of the guanine bases is already close to zero 
and thus the occupation differences are also small. The 
double guanine segment limits the current the molecule 
can support at the adenine energy around — 0.44eV. 

If we artificially enhance the rate Wgg by a factor 100, 



see dashed line on Fig. [9] the conductance is overall in- 
creased, though not by a factor of 100. The increase is 
much stronger around the adenine energy than around 
the guanine energy. Removing the original bottleneck, 
even small occupation differences between the guanines 
lead now to a fairly sized conductance around the ade- 
nine energy. As the conductance at the guanine energy 
increases only slightly despite the 100-fold increased rate, 
the guanine occupation difference at the guanine energy 
is actually reduced almost by the same factor. 

Similar argument hold again for the suppression of the 
adenine related maxima in the differential conductance. 



IV. SUMMARY 

We have presented a diagrammatic real-time approach 
to polaron hopping through molecules coupled to metal- 
lic electrodes, taking into account vibration-mediated 
charge correlations. This technique leads to a hierarchy 
of linear rate equations for the occupation and many- 
particle correlation functions, which is naturally trun- 
cated for a finite size system. Thus, an exact description 
of correlation effects is possible for a given order of the 
perturbation expansion in the hopping parameters. Us- 
ing short DNA molecules as an example, we show that 
including exact correlations lowers the zero bias conduc- 
tance of inhomogeneous DNA sequences when the aver- 
age charge occupation is sufficiently high. For exponen- 
tially small charge occupation, a mean-field description 
of the correlations is adequate. For the I-V characteris- 
tics, the inclusion of exact correlation effects is necessary 
since at the experimentally relevant bias voltage the local 
charge densities are generally non-negligible. We further 
have shown that for short DNA molecules consisting of 
two different types of base pairs, the zero bias and differ- 
ential conductance shows two maxima. Depending on the 
specific sequence, one of these maxima can be suppressed 
due to low hopping rates. 
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Kompctenznetz "Funktionelle Nanostrukturen" . 



APPENDIX A 
Construction of irreducible block diagrams 

Below we will state the rules for the construction and 
evaluation of irreducible block diagrams. The rules for 
pure hopping diagrams, i. e. those containing only ver- 
tices oc tij, were developed by Bottger and Bryksin 12 . 
We extended their theory adding new rules to treat dia- 
grams with tunneling vertices oc t r iv . 

The perturbative expansion can be visualized by the 
construction of diagrams which are equivalent to expres- 
sions in the analytic expansion. The main contribution 
to the diagrams comes from so called irreducible blocks, 
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which, as the name implies, cannot be decomposed into 
more simple diagrams. The main feature of an irreducible 
block diagram is, that it does not diverge, when inte- 
grating over the internal times t$. Irreducible blocks can 
be identified by their property of not allowing free sec- 
tions. A free section is a vertical line drawn between the 
leftmost vertex and the rightmost vertex (except for the 
clamp) that does not cross either a phonon line or an 
external fermion (tunneling) line. 

The rules come in two sets: the first for the construc- 
tion and labeling of possible diagram, the second set for 
the evaluation of a particular diagram. The rules are 
general for all orders of perturbation theory. 

1. Draw the Keldysh time contour as a rectangle 
which is open to the left, corresponding to t — * — oo. 

2. For a diagram of order n we draw on the contour 
n + 1 pair vertices consisting of one open circle 
O (symbolizing a destruction operator) and one 
crossed circle (J) (symbolizing a creation opera- 
tor). All circles belonging to operators acting on 
the molecule are drawn on the inside of the contour, 
whereas circles belonging to electrode operators are 
drawn on the outside of the contour. Therefore, if 
the pair vertex is due to a tunneling process V one 
circle is on the inside and the other one is on the 
outside of the contour. The circles of a hopping 
process are both drawn on the inside of the con- 
tour where the open circle is always 'earlier' along 
the Keldysh contour than the crossed circle. As we 
calculate diagrams to evaluate the density matrix, 
we draw one pair vertex (also called 'clamp'—) at 
the inside of the right vertical line of the Keldysh 
contour, corresponding to time t. The other n ver- 
tices are drawn at n times U on either the upper or 
lower branch of the Keldysh contour (where ti is 
the leftmost, earliest time and t\ is the rightmost, 
latest time). 

3. Each open circles O on the inside of the contour 
has one ingoing fermion line (arrow pointing to the 
vertex) and each crossed circle (J) has one outgoing 
fermion line (arrow pointing away from the vertex) 
which is locally directed along the Keldysh contour. 

4. Complementary circles outside the contour are 
pairwise connected by a fermion line drawn outside 
of the contour. Since this line corresponds to an 
electron propagating in electrode r the connected 
circles have to belong to the same electrode r, oth- 
erwise the diagram contribution is zero. 

5. The clamp is always connected by a fermion line to 
the rightmost vertex (other than the clamp) drawn 
inside of the contour, i. e. the clamp and the right- 
most vertex are associated with the same state. 
If the rightmost vertex is a hopping vertex, the 
fermion line is directed along the contour. If the 
rightmost vertex is a tunneling vertex, the inside 



circle (open or crossed) is connected to the comple- 
mentary circle of the clamp, no matter what the 
direction of the fermion line. 

6. The remaining unconnected inside circles have 
fermion lines going into (coming from) the region 
left of the diagram (t — > — oo) without intersecting 
each other. 

7. Each circle belongs to one specific state for the 
molecule or the electrode. We label the molecule 
states (sites) by latin characters (e.g. m, n, . . .) and 
the electrode states by Greek characters (e.g. v). 
Note that the two circles of a hopping vertex can 
not correspond to the same state (site). Since we 
want to calculate the density matrix pi then both 
circles of the clamp are associated with the state 
(site) I. 

8. Except for the clamp, the circles on the inside of 
the contour must be connected by phonon lines so 
that the diagram has no free section, as defined 
above. One circle can be connected to more than 
one phonon line. All diagrams with different num- 
ber of phonon lines (but still without free sections) 
have to be considered. Only circles belonging to 
the same state (site) can be connected by a phonon 
line. Therefore, the two circles of a hopping vertex 
can not be connected. 

The rules for evaluating a diagram are as follows. 

1. A hopping vertex at time i, is associated with 
a factor ±it ram /C ra /C m e _ ^ e " _c, ™^ ti_ * 2 ' where the 
creation operator (crossed circle) corresponds to 
site label n and the destruction operator to the 
site label m (recall that t% is the leftmost time 
of the diagram). A tunneling vertex is asso- 
ciated with a factor ±it r K m e- l{t »- t ™ ){ - u - t2) or 
±ii r /C n e -I ( e " -e "'(* i-t2 - ) if the creation operator acts 
on the electrode or on the molecule, respectively. 
Vertices on the upper half of the contour have the 
minus sign, vertices on the lower half of the contour 
have the plus sign. The factor 

JC m = exp j-i £ (^j\ 2 N(w a ) + 1) j . 

2. The outside fermion lines of the electrodes r con- 
tribute a factor 1 — or depending whether 
they run in the direction of the contour or against 
it. Here fl is the Fermi function at energy e„ — fi r , 
with the chemical potential [i r . 

3. The fermion lines entering (leaving) the irreducible 
block from (to) the left are labeled from top to 
bottom, with one ingoing and one outgoing line 
belonging to pairwise the same state. The labels 
determine the indices of the irreducible block, e. g. 
W mim2 ;(i). The line attached to the clamp and 
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FIG. 10: 



leaving to the left is associated with a Kronecker 
delta function of the states the line connects. 

A phonon line connecting two circles both associ- 
ated to a state (site) m has a value 

FrkiU -tj) = exp {(A m (ti - tj)} - 1 , with 

. f \ am \ 2 cos(uj a [t + ihf3/2]) 

A m {t) — 



sinh (huj a (3/2) 



where the circle at time U is later on the contour 
than the circle at time tj. The factor £ is deter- 
mined by the type of circles the line connects. If 
the circles are different ( = +1, otherwise £ = — 1. 

4. Multiply with a factor (-1) M+Ar , where M is the 
number of intersections of fermion lines on the out- 
side of the contour (tunneling lines) and N is the 
number of intersections of fermion lines on the in- 
side of the contour. 

5. We integrate over all internal times U (except t\ 
and t 2 ) and sum over all electrode states v and 
all internal molecule states except the states 
associated with the clamp. 

As an example Fig. [TO] shows a second order diagram 
for a hopping process between molecule and metallic elec- 
trode. Following the rules listed above the value of this 
diagram is 

W milJ ={-l){-if\t r \ 2 6 mil Kl 

x^Ue-^-^-^F+iti-h) 



The last line is true in the wide band limit with r L ' R oc 

Pe|f L ' R | 2 . 

APPENDIX B 

Determination of the correlation function 

The correlation function with which the irreducible 
block is convoluted can be easily constructed from the 
diagrams. Lines leaving the block correspond to creation 



operators, lines entering the block correspond to destruc- 
tion operators. The order of operators from left to right 
in the correlation functions corresponds to the order of 
the lines leaving/entering the diagram from bottom to 
top. 

For example, let us consider some irreducible block 
W miTO2 m 3 ;(t) which has the following order of the ter- 
minals at the left of the diagram from bottom to top: 
TO3 outgoing, m 2 ingoing, 7713 ingoing, mi outgoing, m 2 
outgoing, and mi ingoing. Thus this rate will be convo- 
luted with the third order correlation function 



U £ d (0)c4 3 a m2 a m3 aj^ a) m2 a mi U ^ {0)J, = 

Pmim 3 — Pm 1 m 2 m 3 

APPENDIX C 
Vibrational operator products 



In the perturbation expansion of the single particle 
density matrix pi(t) = ^aj(t)ai(t)^ . to order n in the 

perturbative Hamiltonian H' (Eq. [5|) , one obtains up to 
n vibrational operators (equal number of x and X ) at dif- 
ferent times which act upon the same vibrational states. 



Xk{ti)x k (ti)Xk{h) ■ ■ ■ X k (tn) 



exp 



Afcn 



x exp { {(uAkih - t 2 )} + T c {(i 3 A k {ti - t 3 )} 



+T C {{3,nMh-t n )} + 



1 when Xk(ti)xl(tj) or xUWXkiU), 
1 when Xk(U)xk(tj) or X f k (U)xl(tj) 

The expression Tc in Tc {(i 2 A k (ti — ^2)} ensures, that 
ti is later on the contour than t 2 and A k (ti — t 2 ) is given 
by 




A k {t 1 -t 2 ) 



E 



Xka 
U ka 



cos (u ka [ti - t 2 + ih/3/2]) 
sinh (huj ka [3/2) 



For a correlator with n operators Xk and x\ acting on 
the same state one gets N — C")^" 1 ) different terms 
A k (ti — tj) in the exponential function. This is due to 
the various operator commutations involved in deriving 
the above expression. 
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